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ABSTRACT 

An  evaluation  is  performed  of  a  radiation  model  for  the 
Naval  Postgraduate  School  primitive-equation  numerical  wea- 
ther prediction  system.   The  model  employs  empirical  expres- 
sions for  atmospheric  absorptivity,  scattering-reflectivity, 
cloud-reflectivity  and  earth-surface  reflectivity  to  compute 
solar  insolation   absorbed  at  earth  and  in  the  key  atmos- 
pheric layers.   The  terrestrial  cooling  effect  at  earth  and 
in  these  same  key  atmospheric  layers  is  formulated  using  re- 
cent empirically-derived  emissivities  for  the  effects  of 
both  water  vapor  and  C0„ . 

Mean  seasonal  atmospheric  soundings  for  the  Northern 
Hemisphere  are  utilized  for  testing  the  model.  In  addition, 
application  of  atmospheric  boundary-layer  modeling  permits 
determination  of  the  surface-layer  turbulent  transports  (in 
the  vertical)  of  sensible  and  latent  heat  at  the  earth's 
surface.   To  evaluate  the  validity  of  the  radiation  model, 
heat  budgets  are  compiled  for  the  earth-atmosphere  system, 
the  atmospheric  column  and  the  earth's  surface. 
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I.   INTRODUCTION 


The  present  study  was  to  evaluate  a  radiation  package 
for  the  Naval  Postgraduate  School  primitive  equation  pre- 
diction model.   The  radiation  package  is  primarily  an  out- 
growth of  the  model  used  in  the  UCLA  general  circulation 
model  (Gates  et  al.  1971)  with  modifications  suggested 
from  the  NCAR  model  (Oliger  et  al.  1970)  and  from  the  Geo- 
physical Fluid  Dynamical  Laboratory  (Manabe  and  Strickler, 
1964)  . 

The  model  consists  of  five  a-levels  (a  =  P/tt  where  77 
is  the  surface  pressure)  for  the  purposes  of  data  input 
as  shown  in  Figure  1.   Cloud  cover  was  treated  by  the  use 
of  a  single  cloud  layer  extending  from  a  =  0.8  to  0.4. 
In  the  computations  for  both  pressure-scaled  water-vapor 
and  CO-  absorber  masses,  a  more  realistic  method  was  in- 
troduced to  extrapolate  accurate  water-vapor  values  above  . 
levels  normally  available  from  radiosonde  data  (P  ~  300  mb) . 
A  system  for  the  integrations  of  net  terrestrial  flux  was 
introduced  by  using  empirically  derived  absorptivit ies 
(Sasamori,  1968)  based  on  the  radiation  chart  of  Yamamoto 
(1952),  for  water-vapor  and  C0_ .   Evaporation  and  sensible 
heat  fluxes  were  derived  from  a  formulation  for  surface 
layer  heat  transfer  by  Langlois  and  Kwok  (1969). 

Northern  hemisphere  climatological  atmospheric  data 
from  London  (1957)  were  utilized  to  evaluate  the  model. 


Pressure,  temperature,  and  relative  humidity  data  were 
divided  into  the  four  seasons  and  available  over  ten  de- 
gree latitude  bands  from  the  equator  poleward. 

The  development  of  the  terrestrial  fluxes  is  given 
in  Section  III  and  the  solar  radiation  model  is  given  in 
Section  IV.   The  results  of  the  computed  long  and  short 
wave  fluxes  were  then  compared  to  those  of  Budyko  (1955), 
London  (1957)  and  Katayama  (1967). 

Finally,  to  evaluate  the  adaptability  of  the  radiative 
package  to  a  general  circulation  model,  three  heat  bud- 
gets were  attempted.   The  budgets  were  for  the  earth-at- 
mosphere system,  the  atmospheric  column  and  the  surface. 
Since  heat  transfer  calculations  due  to  evaporation  and 
sensible  heat  are  more  easily  made  over  ocean  surfaces, 
London's  atmospheric  data  was  fitted  to  climatological 
ocean  surface  temperature  and  wind  data  from  the  United 
States  Weather  Bureau's  Atlas  of  Climatic  Charts  of  the 
Oceans  (1938).   Surface  layer  evaporation  and  sensible 
heat  flux  computations  were  not  conducted  poleward  of  60N 
because  of  the  possibility  of  ice  cover.   Rather,  com- 
bined evaporation  and  sensible  heat  rates  were  taken  from 
Vowinckel  and  Taylor  (1964). 
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Figure  1.   Representation  of  the  five-layer  sounding 
model  used  for  radiative  transfer  calculations-.   Tem- 
perature, mixing  ratio,  pressure-scaled  water-vapor 
(Au)  and  C0?  (Ac)  absorber  mass  increments  are  shown 
a.t  the  appropriate  G-levels  and  k-levels  (k  =  10a)  . 
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II.   DATA  PREPARATION 

A.   TEMPERATURE  AND  MOISTURE  TERMS 

Northern  hemisphere  vertical  and  latitudinal  distribu- 
tions of  pressure,  temperature,  and  relative  humidity  were 
taken  from  the  climatological  soundings  of  London  (1957). 
The  data  was  compiled  in  four  mean  seasonal  tables.   The 
latitudinal  variation  consisted  of  10  degree  latitude  bands 
measured  from  the  equator  northward.   In  the  vertical,  data 
was  distributed  at  1  km  intervals  from  the  surface  to  the 
base  of  the  stratosphere. 

With  actual  radiosonde  data,  continuous  information 
would  be  available  above  the  tropopause  to  at  least  the 
0.1  sigma  level.   The  a    =  0.1  level  was  the  highest  level 
at  which  the  water-vapor  mixing  ratio  was  considered  to  be 
recoverable  from  the  fields  of  q.   London's  radiosonde  data 
was  continuous  only  to  the  tropopause.   Since  the  tropo- 
pause varies  in  height  latitudinally  as  well  as  seasonally, 
adjustment  of  the  vertical  distribution  of  T  and  q  had  to 
be  made  to  insure  that  the  data  would  be  available  to  a  = 
0.1  level.   The  simplest  solution  for  the  stratosphere  was 
to  assume  isothermal  conditions  with  constant  relative 
humidity  to  the  a  =  0.1  level  for  any  vertical  sounding 
that  ended  below  a  =  0.1. 

It  was  necessary  to  interpolate  from  the  distribution 
at  successive  km  levels  in  the  vertical  (as  listed  in 
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London,  1957)  to  the  specific  sigma  levels  of  Figure  1. 
Temperatures  were  carried  at  the  1.0,  0.8,  0.6,  0.4,  0.2 
and  0.1  a-levels.   Since  surface  data  was  available  in 
London's  soundings,  Tlf.  was  set  equal  to  the  listed  sur- 
face temperature  rather  than  requiring  an  extrapolation 
of  the  type 

T1Q  =  0.5  (3T9  -  T?)  (1) 

Temperatures  at  other  sigma  surfaces,  say  at  level  k,  were 
interpolated  using  a  three-point  Lagrangian  formula  of  the 
form 

T(k)  =  (T2)(Pk-P1)(Pk-P0)/(P2-P1)(P2-P()) 

+  (T1)(Pk-P2)(Pk-P0)/(P1-P2)(P1-PQ) 

+  (V(Pk-Pl)(Pk-P2)/(P0-P2>(P0-V    (2) 

where  it  is  assumed  that  the  data  level  (P  ,T  )  lies  above 

the  a-level  in  question  and  (P_,T_)  lies  below  that  level. 

Mixing  ratio  values  were  calculated  with  London's  data 

using 

q  -  RH  (0.622e  ) /P  (3) 

s 

The  q-values  were  then  interpolated  to  a, -levels  using  the 

Lagrangian  formulation  of  (2)  applied  to  q's  at  known 

levels.   The  q,  -values  were  centered  at  the  odd  a-levels 
k 

of  0.9,  0.7,  0.5  and  0.3.   These  values  are  then  considered 
to  be  the  layer  means  between  the  even  a-levels. 

Accurate  humidity  information  is  in  general  not  avail- 
able above  a  =  0.3,  especially  in  the  stratosphere.   To 
find  a  mean  q-value  representative  of  the  layer  (0,2)  a 
different  formulation  was  used. 
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Smith  (1966)  has  shown  a  strong  correlation  between 
q(P)  and  the  surface  value  q ..  „  through  the  ex  t  rapola  t  ive 


formula 


lill  =  (JL_) 


x 


(4) 


^10       10 
London's  soundings  were  tested  to  obtain  a  best-fit  be^ 
tween  In  q  and  In  P  of  the  form 


sin  .  (jl.}x 

q3       P3 


(5) 


which  is  similar  in  form  to  that  suggested  by  Smith.   The 
values  of  q(P)  resulting  from  (5)  for  a  ■  0.9,  0.7,  0.5 
and  0.3  were  used  to  solve  for  the  profile  parameter  X. 
This  parameter  varied  with  latitude  and  season  as  shown 
in  Table  I.   The  relationship  (5),  with  appropriate 
X-values,  was  then  used  to  extrapolate  q  at  a  =  0.1.   The 
resulting    q  -values  were  then  considered  the  mean 
q-values  for  each  layer  (0,2). 

The  procedure  just  described  to  derive  q     was  re- 
garded as  justifiable  in  view  of  the  following  reasons: 

(i)   The  residuals  [q(est)  -  q(obs)],  were  small  for 
all  k  ^  3  ,  and 

(ii)  The  extrapolated   q    corresponds  to  the  proper 
order  of  magnitude  as  listed  in  model  atmospheres  (Valley, 
1965)  at  pressures  P  ~  100  mb . 
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TABLE  I.   Seasonal  and  latitudinal  values  of  the  profile 
parameter  X  of  Equation  (5). 


Latitude  Winter  Spring  Summer  Fall 

0-10  3.1438  3.0178  3.2885  3.0202 

10-20  3.1182  3.3024  3.4381  3.2555 

20-30  3.1205  3.3978  3.4272  3.2374 

30-40  2.6437  3.3838  3.3563  3.0873 

40-50  3.0002  3.5178  3.2480  3.1259 

50-60  2.7210  3.2382  3.2031  3.1846 

60-70  2.3267  2.7478  3.0684  2.8925 

70-80  1.9671  2.1953  2.7886  2.4617 

80-90  1.1743  1.7831  2.8434  2.1262 
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B.   CLOUD  MODELING 

In  both  the  solar  and  terrestrial  radiative  transfer 
processes,  it  is  essential  to  know  the  fractional  cloud 
cover  (CL)  at  each  grid  point.   Clouds  act  as  almost  perfect 
blackbody  radiators  for  infrared  radiation,  and  as  efficent 
reflectors  for  solar  radiation. 

There  is  a  difficulty  in  parameterizing  different 
cloud  types  and  heights.   In  this  model,  as  in  the  UCLA 
general  circulation  model,  only  a  single  cloud  layer  was 
parameterized  for  the  purposes  of  inclusion  in  the  vertical 
radiative  transfer  model  at  gridpoints.    As  is  done  in 
most  other  numerical  prediction  schemes  involving  cloud 
amount,  CL ,  use  is  made  of  Smagor insky ' s  (1960)  algorithm 
(applicable  at  approximately  700  mb) 


CL  =  2.0  (e/ee)  -  0.7 


(6) 


where  0  <^  CL  <^  1.0.   When  a  cloud  amount  CL  >  0  existed, 
it  was  assumed  to  extend  between  a  =  0.8  to  a  =0.4. 

Comparison  of  the  calculated  values  of  CL  computed 
from  (6)  and  those  tabulated  from  seasonal  climatological 
estimates  by  Telegadas  and  London  (1954)  and  Seinde  (1954) 
showed  a  marked  difference  in  the  mid-latitudes  as  seen 
by  comparison  of  Tables  II  and  III.   Since  other  climato- 
logical estimates  of  cloud  cover,  namely  Sellers  (1965), 
compared  best  with  the  data  tabulated  in  Table  II, climato- 
logy was  used  in  this  model  rather  than  the  results  of 
Table  III. 
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TABLE  ii.  The  seasonal  distribution  of  total  cloudiness 
(percent  cloud  cover)  tabulated  by  London  (1957). 


Lat  N. 

Winter 

Spring 

Summe  r 

Fall 

0-10 

47 

51 

54 

53 

10-20 

36 

42 

49 

48 

20-30 

38 

42 

42 

41 

30-40 

50 

52 

41 

46 

40-50 

59 

59 

55 

56 

50-60 

63 

62 

63 

66 

60-70 

58 

60 

66 

70 

70-80 

47 

59 

69 

70 

80-90 

40 

55 

64 

60 
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TABLE  III. The  seasonal  distribution  of  total  cloudiness 
(percent  of  sky  cover)  calculated  by  the  Smagor insky 
algorithm. 


Lat  N 

Winter 

Spring 

Summer 

Fall 

0-10 

30 

40 

56 

46 

10-20 

12 

20 

.  48 

36 

20-30 

10 

10 

30 

20 

30-40 

18 

14 

22 

16 

40-50 

34 

26 

30 

26 

50-60 

42 

34 

42 

42 

60-70 

46 

40 

50 

46 

70-80 

46 

38 

50 

48 

80-90 

34 

28 

50 

40 
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Even  though  London's  cloud-cover  climatology  was 
accepted  as  valid  here,  it  should  be  realized  that  London's 
results  included  overlapping  clouds  at  multiple  levels, 
which  was  impossible  with  the  parameterization  using  (6). 
An  improved  cloud  model  permitting  cloud  existence  at 
two  or  more  levels  in  the  vertical  should  be  attempted 
in  future  radiative  heating  systems.   At  the  present  time, 
the  added  complexity  of  the  radiative  physics  is  considered 
too  costly  in  computer  time  for  inclusion  in  the  present 
model.   However,  it  should  be  noted  that  Arakawa  et  al. 
(1972)  have  already  developed  the  concept  of  a  three-layer 
radiative  cloud  model  for  a  1972  version  of  the  UCLA 
general  circulation  model. 

Another  reason  for  the  choice  of  the  climatological 
CL  must  be  considered.   With  the  use  of  climatology, 
transient  synoptic  scale  perturbations  and  disturbances 
in  the  soundings  are  smoothed  or  eliminated.   Therefore, 
the  day-to-day  variations  of  relative  humidity  tend  to  be 
smoothed  out  during  the  season.   The  resulting  smoothed 
relative  humidity  distributions  thus  tended  to  be  less 
than  their  normal  peak-values  and  give  too  small  a  magni- 
tude of  CL  when  employing  the  Smagorinsky  algorithm. 

C.   ABSORBER-MASS  DEFINITIONS:   PRESSURE  SCALED 

Since  mixing  ratio  values  have  now  been  computed  at 
the  five  odd  O-levels  (Figure  1),  it  is  possible  to  deter- 
mine the  pressure-scaled  water-vapor  absorber  masses, 
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Au,  ,  associated  with  the  five  corresponding  layers.   This 
is  done  in  the  model  in  such  a  way  that  the  additive 
scheme 

u(6,10)  =  u(8,10)  +  u(8,6)  (7) 

holds.   The  layer  water  vapor  absorber  mass,  u(2k,  2k-2) , 
has  been  centered  at  odd  levels  (2k-l) .   The  scaled 
absorber  mass  is  calculated  using 


Au 


q2k-l  AP  ,P2k-lN0.72 


2k-l 


( 


(8) 


8        '0 
which,  as  has  been  noted,  is  the  absorber  mass  in  the 

layer  (2k,  2k-2) .   AP  is  a  constant  in  the  analysis  scheme 

equal  to  0.2ir.   The  mixing  ratio  value,  q9v_-,»  is  the 

mean  centered  value  in  the  layer  (2k,  2k-2) .   The  pressure 

0.72 
ratio  factor  (P„,_1/P0)      is  the  scaling  factor  for 

collisonal  line-broadening  after  Moller  and  Raschke  (196A) 

This  latter  formulation  disposes  of  the  necessity  of 

temperature  scaling  in  the  usual  Lorentz  line-broadening 

formula  (c.f.,  Danard,  1969),  where  the  pressure-ratio 

scaling  is  linear  and  a  temperature-ratio  scaling  factor 

is  indicated  by  (T/TQ)~  '   times  the  N.T.P.  mass  within 

a  layer.   The  subscript  "0"  denotes  an  N.T.P.  condition 

where  Pn  =  1013.25  mb ,  Tn  =  273.16°K,  and  pn  is  the  N.T.P. 


0 


value  of  density 


0 


0 


If     (8)     is    integrated    from    P2       =    (0.2k)    PQ    to    P    =    ir , 

the    result    is 

P  n=k> 1 

t    /t>    n!.72/    r0         "    £-      -  r    1.72  1.72,       ,Qv 

U2k    =    (7r/V  (l77Tg)       Z,       «2n-l     [a2n         "    °2n-2]        (9) 

n  =  j 
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where  by  definition 

P    P 

ak  "  F2k/7T  "  PQ    7T  (10) 

For  CO2,  Moller  and  Raschke  (1964)  proposed  pressure- 
scaling  similar  to  that  of  water-vapor  but  with  a  pressure- 
ratio  exponent  of  0.65.   For  simplicity,  a  single 
pressure-scaling  exponent  of  0.72  was  utilized  for  both 

cases.   The  form  of  the  layer-reduced  absorber  mass  of  C0_ 

2 
in  N.T.P.  cm/cm   in  the  layer  (2k,  2k-2)  is 


p 
Ac2k_1  =  3.14  x  10-VAP/gp0V-P-=i)0-72         (11) 


p 
0 

The  process  of  integration  from  P_,  to  P~  yields  the 

formulation 

-4        P    1.72 
_  3.14  x  10     r      2k   *   ,  ,10v 

c2k  =  1772 [     V      ]  (   } 

where  H  =  7.995  x  10  cm  is  the  height  of  the  homogeneous 
atmosphere  in  N.T.P.  cms.   Adjusting  (12)  to  a  O-coordinate 
system  with  P  =  ij  at  the  surface  yields 

3.14    x    10"4    u,Tr    N1.72M       n      1.72.  .... 

c2k  =  ITTi h(p"0)  [1"  a2k         ]  (13) 

Table  IVshows  the  seasonal  distribution  of  the  total 
pressure-scaled  water-vapor  absorber  masses,  u(0,10).   As 
expected,  the  summer  season  contains  the  largest  absorber 
mass  with  the  least  in  winter.   Spring  and  fall  values  are 
intermediate  between  the  maximum  and  minimum  values.   Com- 
putations for  the  C0_  absorber-masses  are  dependent  upon 
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TABLE  IV.  Seasonal  distributions  of  total  pressure-scaled 
water-vapor  masses,  u(0,10),  in  gm  cm^  . 


Latitude 

Spring 

Summer 

Fall 

Winter 

0-10 

3.878 

4.287 

4.023 

3.624 

10-20 

3.271 

4.316 

3.826 

2.749 

20-30 

2.417 

3.601 

3.055 

1.993 

30-40 

1.607 

2.897 

2.133 

1.279 

40-50 

1.099 

2.292 

1.426 

0.747 

50-60 

0.703 

1.793 

0.932 

0.449 

60-70 

0.437 

1.428 

0.598 

0.256 

70-80 

0.289 

1.045 

0.471 

0.163 

80-90 

0.192 

0.901 

0.319 

0.105 
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the  total  mass  of  air  in  a  unit  column,  and  hence  on  the 
surface  pressure.   Since  the  variation  in  surface  pressure 
averaged  across  all  latitude  bands  was  less  than  10  mb , 

the  CO   absorber-masses  were  nearly  constant  [c(0,10)  = 

2 
125  cm/ cm  ] . 
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III.   TERRESTRIAL  RADIATIVE  TRANSFER 

A.   THEORETICAL  AND  EMPIRICAL  BACKGROUND 

Virtually  all  of  the  radiation  packages  in  present  use 
with  large  scale  circulation  models  in  the  United  States 
[cf.,  Manabe  and  Strickler  (1964),  Oliger  et  al.  (1970), 
Arakawa  et  al.  (1969)]  make  use  of  Yamamoto's  (1952)  water- 
vapor  flux  emissivities  as  a  basis  for  their  long-wave 
radiative  flux  calculations.   The  water-vapor  flux  emissi- 
vity  is  normally  thought  of  as  a  function  of  both  water- 
vapor  scaled  mass  and  temperature.   Yet,  in  the  normal  range 
of  temperatures  (T  >_  220K) ,  the  emissivity  is  virtually 
independent  of  temperature.   At  temperatures  T  <  220K,  the 
wave-averaged  flux  emissivity  is  defined  after  Yamamoto 
(1952)  by 

'  [1"TFV(1Vou)]  7T-(T=220K>  dV 

e(u,T)  =  -2^ (14) 

r   dB 

Here  the  numerator  has  its  temperature  dependent  terms 

3 
fixed  at  220K, while  the  denominator  equals  4St(T  )  even 

when  K220K.   The  temperature  affects  the  population  of 

different  rotational  energy  levels  within  the  so-called 

"rotational"  band.   At  low  temperatures,  Planck's  function 

shifts  maximal  radiant  energy  to  the  longer  wavelengths 

in  the  rotational  band,  according  to  Wien's  displacement 
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law  of 

X  =  2897/T  (15) 

max 

where  A     is  in  microns.   Thus  with  T  <  2  2 OK,  the  highly 

max  °       J 

absorbed  rotational  lines  are  the  only  contributors  to 
water-vapor  absorptivity.   This  leads  to  an  increase  of 
absorption  near  the  center  of  the  rotational  band  and  thus 
to  an  increase  in  mean  absorptivity  with  decreasing  tem- 
perature.  This  phenomenon  is  depicted  graphically  in 
Figure  2  after  Sasamori  (1968). 

Yamamoto's  radiation  chart  has  been  the  basis  of  most 
recent  radiation  calculations.   Sasamori  (1968,  1970)  used 
Yamamoto's  radiation  theory  to  write  empirical  formulas 
for  all  of  the  emissivities  needed  for  long-wave  radiative 
flux  calculations.   These  empirical  formulas  have  been 
shown  by  Sasamori  to  yield  adequate  accuracy  and  were 
adapted  for  this  numerical  study.   The  temperature  indepen- 
dent emissivity  for  water-vapor,  from  Sasamori,  is  given  by 

e(u)  =  0.24  log1()(u  +  0.010)  +  0.622  (16) 

This  formulation  was  used  for  T  •>  220K.   For  T  <  220K, 
Sasamori  gives  the  temperature-dependent  part  of  the  emis- 
sivity function  by  water  vapor  as 

e(u>I)=(8.34T0-3531o«10u-0-44)u-0-34551o8iou--705    (17) 

When  absorption  by  C0»  and  water-vapor  are  considered 
jointly,  the  product  law  of  transmissivit ies  must  be 
employed.   In  doing  this,  the  transmiss ivity  of  the  two 
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e(u) 


K> 


u(gm  cm   ) 


Figure  2.   Temperature  dependent  part  of  the  emissivity 
with  varied  amounts  of  scaled  water-vapor. 
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absorbers  acting  jointly  may  be  expressed  as 
TF(u,c)  =  Tf(u)tf(c) 

The  flux  integral  as  formulated  by  Yamamoto  is 


(18) 


u(T) 
F  -   /       e(u,T)  dB 

U=0'TRef 


(19) 


where  u(T)  is  the  water-vapor  sounding  mass  measured  away 
from  the  reference  level.   Equation  (19)  then  becomes  with 
the  inclusion  of  C0„ 


u,c(T) 
F  =    / 

u=0,Q=0 


[l-TF(u)TF(c)]dB 


(20) 


or 


u.c 


F  »    f  [e(u,T)  +  Ae]  dB 

u=0,c=0 


(21) 


The  Ae  term  represents  the  effective  emission  due  to  CO 
not  already  counted  in  the  water-vapor  flux  calculations. 
The  added  emission  is  expressed  on  the  Yamamoto  chart  by 
an  Increase  of  emissivity  at  all  points  of  the  radiative 
sounding  (except  at  the  reference  level  points).   The  Ae 
term  includes  only  the  15um  band  emitted  flux  of  C0„  which 
has  not  already  been  counted  in  emission  by  overlap  with 
the  rotational  band  of  water  vapor.   Sasamori  considered 
this  C0„  band  emissivity  to  be  temperature  independent 
in  developing  the  empirical  formulation 

Ae  =  0.07262[l-0.62556(u+0.286)0,26] [logc+1.064]   (22) 
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In  the  following  flux  computations,  the  emissivity,  e , 

will  include  both  the  effects  of  C0?  given  by  (22)  and  that 

of  water  vapor  given  in  (16)  and  (17). 

B.   TECHNIQUES  OF  NET  FLUX  COMPUTATIONS 

Schematic  representations  of  the  Yamamoto  radiation 
chart  were  utilized  for  considering  the  integration  of  the 
net  long  wave  flux  calculations.   In  Figure  3,  the  hatched 
area  is  a  representation  of  downward  flux  at  the  surface. 
The  area  may  be  computed  using  the  trapezoidal  summation 
rule  as  follows: 

Ff(XO)  =  0.5{[e(8-10)(B10-B8)]+[e(8-10)+e(6-10)] 
(Bg-B6)+[e(6-10)+e(4-10)] (B6-B4)+[ e (4-10) 
+  e(2-10)](B4-B2)+[e(l-10)+  e(2-10)] 

B(T  ) 
(B2-B1)+  2  /    X   E[(0-10),T]db} 

B=0  .    (23) 

Here  B  =  St  *  T   is  the  S tef an-Bolt zmann  blackbody  flux 
arid  e(8-10),  for  example,  represents  the  emissivity  com- 
puted from  Sasamori's  empirical  equations  using  the  water 
vapor  and  CO   present  in  the  layer  (8-10) .   The  term 
£  =  e(u,T)  +  A£  is  solved  by  e(u,T)given  by(17)and  Ae  by 
(22).   The  subscripts  on  the  S t ef an-Bolt zmann  fluxes  (i.e., 
B.Q-  Bg)  are  identifers  for  the  layer  contribution  from 
which  the  flux-emission  takes  place. 
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Figure  3. 
level  10. 


Schematic  area-depiction  of  net  flux  at 
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B(T  ) 
The  final  integral  of  (23),   /      e  dB  , 

B=0 

may  be  evaluated  by  the  integrated  result 


St    *    t/    {E(uo,  T    )(___ 


g10    uo+    3.56 
0.26 


)    + 


0.07262[l-0.62556(u    +    0 . 286)  ] [log. ftc    +    1.064]} 

o  10  o 

(24) 

In  (24),  it  is  understood  that  the  integration  spans 

the  temperature  range  (0,T.)  along  the  final  isopleths 

u   and  c   (e.g.,  Figure  3).   The  net  flux  at  level  10  is 
o       o 


given  by 


Fio  "  Bio  -  F*(10> 


(25) 


where  the  asterisk  superscript  denotes  F(up)  -  F(down)  at 
the  level . 

The  net  flux  at  level  6  is  represented  schematically 
by  B..-  minus  the  hatched  area  in  Figure  4.   The  numerical 
modeling  follows  (23)  where 

F6*  =  B1Q  -  0.5  U(6-8)(B8-  Bfi)  +  [e(6-8)+  e(6-10)] 

(B1Q-  B8)  +  £(4-6)(B6-  B4)  +  [e(4-6)+  e(2-6)] 

(B4-  B2)  +  [£(2-6)  +  e(l-6)](B2-  B1)  + 


B(T  ) 

2   /      e[(0-6),  T]  dB} 
B  =  0 


(26) 
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Figure    4. 
level    6 . 


Schematic    area-depiction    of    net    flux    at 
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Figure  5  is  the  schematic  for  net  flux  at  level  2. 
Following  the  previous  computations 

F2*  =  B10  "  °*5  {e(2~4)(B4-B2>+[£<2-4)+  £(2-6)] 
(B6-  B4)+[e(2-6)+  e(2-8)](Bg-  Bfi)  + 
[e(2-8)+  e(2-10)](B10-  Bg)  +  e(l-2)(B2~  B.^ 


B(T  ) 
+    2       f       X       e     [  (0-2)  ,  T]  dB} 
B  =  0 


(27) 


The  foregoing  net  fluxes  F    ,  F,  ,  F_   are  clear  sky 
calculations  only.   A  modification  due  to  cloudy  skies  is 
made.   The  cloud-covered  fraction  is  given  by  CL  and  the 
clear  sky  fraction  is  (1-CL) .   The  cloud-covered  fraction 
is  considered  an  overcast  which  extends  vertically  from 
level  8  to  level  4.   Figure  6  shows  a  schematic  flux  dia- 
gram, after  Yamamoto  (1952),  of  net  fluxes  at  levels  10 
and  2,  when  only  the  overcast  sky  fraction  is  considered. 
The  horizontally  hatched  area  is  now  the  net  flux  at  level 
10.   The  net  flux  at  level  6  is  zero  since  a  layer  of  ap- 
proximately 50  meters  in  depth  with  average  cloud  droplet 
content  may  be  regarded  as  a  blackbody  (Brunt,  1939)  at 
the  mean  temperature  of  the  cloud  layer.   The  net  flux  at 
level  2  is  the  slant  hatched  area  in  the  diagram. 

The  net  flux  at  level  10,  with  overcast  skies  is  given 
by 

F10*  =  (B10  "  B8)(1  "  °-5£(8-10))  (28> 
Normally  (28)  represents  a  cooling  effect  at  the  earth's 
surface.   At  the  same  time,  the  net  flux  at  level  8,  or  the 
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Schematic  area-depiction  of  net  flux  at 
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Figure  6.   Schematic  area-depiction  of  net  flux  at 
reference  levels  10  and  2  when  an  overcast  clou.d 
extends  between  levels  8  and  4. 
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cloud  base,  is  an  equal  magnitude  (Fig.  6)  so  that  the 

layer  (8-10)  has  zero  flux  divergence.   However,  the  cloud 

base  receives  a  net  flux  warming  equal  in  magnitude  to  F  _ 

The  net  flux  at  level  2  for  the  overcast  case  is  given 

by 

F2*  =  B^  -  0.5{e(2-4)  (B4-  B2)+  e(l-2)(B2>  B^ 

B(T  ) 
+  2   /  X       e[(0-2),  T]  dB} 

B=0  (29) 


This  net  flux,  F_  ,  is  much  smaller  than  that  for  clear 
skies  because  of  the  lead  term  B,  on  the  right  side  of 
(29)  as  contrasted  with  B  _  in  (27). 

The  net  fluxes,  computed  at  the  three  reference  levels, 
must  now  be  cloud-weighted  to  yield  final  results  for  net 
long-wave  flux  transfer  across  the  interface  by 

F10*  =  tF*o<28>]  CL  +  [F*0(25)](1-CL)        (30) 


* 


* 


F6    =  [F£(26)]  (1  -  CL)  -       (31) 

V   =  tF2(29)]  CL  +  [F2(27)](1-CL)  (32) 

The  temperature-dependent  flux  integral  denoted  by 

■  /  £  (u  ,  Tn )  dB      was  evaluated  for  an  entire  range 
o    .JL 

of  final  Tn  and  u   values  which 
1       o 

could  occur  in  the  atmosphere.   These  results  are  listed 
in  Appendix  Table  A-l. 

C.   DIVERGENCE  OF  NET  FLUX  AND  RESULTANT  COOLING  RATES 

The  resultant,  cloud  weighted,  seasonal  distributions 
of  net  flux  at  the  three  reference  levels  of  a  =  0.2,  0.6, 


34 


and  1.0  are  listed  in  Tables  V(a-d).   The  24-hour  tempera- 
ture changes  between  the  reference  levels  due  to  terrestrial 
radiation  alone  are  also  listed.   The  mean  radiative  tem- 
perature change  centered  at  the  midpoint  of  the  column  was 
computed  by  the  first  law  of  thermodynamics  so  that  the  cool- 


ing rate  for  the  layer  (6,10)  is 


3t    u6 


and  for  the  layer  (2,6)  is 


_,   *  ._  .   1000 

F    )/0.4lT  c 

10  g    P 


3T    ,w  *    v    *  ,n    ,       1000 
¥F-  (F2   -  F6  )/0.4ir  — j-  cp  , 


(33) 


(34) 


where  7T  is  surface  pressure  in  mb  . 

To  compare  the  accuracy  of  the  long-wave  calculations, 
the  meridionally  averaged  mean  values  for  the  fluxes  at  the 
levels  a  =  0.2  and  1.0  were  computed.   The  result  was  repre- 
sented by  the  percentage  of  the  meridionally  averaged  inso- 
lation at  the  top  of  the  troposphere.   Katayama's  63  percent 
(1967),  and  London's  62  percent  (1957)  for  net  flux  at  200mb 
are  in  good  agreement  with  the  64  percent  from  this  study. 
The  net  flux  at  the  surface  was  18  percent  for  all  three 
studies.   Converting  these  spatially  averaged  results  to 
annual  means  (ly  day   ),  we  get  a  net  outgoing  radiation  at 
P  ~  200  mb  of  430  from  the  Katayama-London  means  compared 
to  our  450.   At  the  surface,  the  net  outgoing  long-wave  radia- 
tion is  120  (ly  day   )  compared  to  our  131. 

In  the  net  flux  calculations,  a  comparison  was  made  be- 
tween three  different  atmospheric  situations  in  order  to 
evaluate  the  contribution  of  cloudiness  to  the  total  outgoing 
flux  from  the  earth-atmosphere  system.   One  test  considered 
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only  clear  skies,  another  had  average  cloudiness  and  the 
third  tested  totally  the  overcast  model.   Results  were  com- 
piled in  mean  annual  latitudinal  distributions  of  long- 
wave net  flux  in  Tables  Vl(a-c). 

Since  temperatures  normally  decrease  with  height  in 
the  atmosphere,  a  positive  net  flux  occurred  at  all  levels. 
Strong  surface  temperature  inversions  can  lead  to  negative 
F  n  ,  but  these  conditions  were  not  present  in  this  study. 
Also,  in  an  atmosphere  where  the  effective  optical  depth 
decreased  sharply  with  height,  the  net  radiative  flux  in- 
creases sharply  with  height  at  the  levels  of  the  decreased 
water-vapor  mass.   This  happens  primarily  because  most  of 
the  emitted  flux  from  the  lower  layer  passes  upwards   un- 
attenuated  by  the  dry  upper  layers. 

The  net  flux  through  the  tropopause  is  larger  at  all 
latitudes  for  clear  skies  than  for  conditions  of  average 
or  overcast  cloudiness.   This  occurs  because  the  net  outgo- 
ing radiation  comes  from  a  lower,  and  therefore  warmer 
layer  in  the  case  of  clear  skies.   In  the  case  of  overcast 
skies,  the  radiating  layer  is  the  cloud  top  at  O    =  0.4 
Rather  than  a  =  1.0  for  clear  skies.   The  tropopause  net 
flux  with  average  cloudiness  is  a  weighted  average  between 
the  overcast  and  clear-sky  cases  of  Table  VI. 

The  largest  variability  of  net  flux  magnitudes  with  cloud 
cover  occurred  at  the  O    =  0.6  level.   The  minimum  value  of 
zero  is  observed  during  total  cloudiness  since  this  level  lies 
within  the  cloud  layer  which  was  treated  as  a  blackbody.  Maxi- 
mum values  occurred  of  course  in  conditions  of  clear  skies. 
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IV,   SOLAR  RADIATION  TREATMENT 


A.   EFFECTIVE  SOLAR  INSOLATION 

The  solar  radiation  treatment  is  not  as  straightforward. 
Following  Joseph  (1966),  a  solar  constant  of  2.00  ly  min 
was  assumed.   A  4%  attenuation  by  oxygen  and  ozone  above 
the  tropopause  was  assumed  which  gave  an  effective  solar 
constant  of  1.92  ly  min   .   The  effective  tropopause  height 
was  assumed  to  be  at  a  =  0.2. 

The  effective  insolation  at  the  tropopause  was  computed 

from 

-2 


F  =  S(r/r  )     D  cos  6  (35) 

m 

where 

5  =  effective  solar  constant 

r/r   =  radius  vector  of  earth  relative  to  the  sun 
m 

6  =  zenith  angle  of  sun  relative  to  observing  point 


cos  6=  daily  average  of  cos  9 

D     =  fractional  length  of  daytime 

The  (r/r  )  is  listed  as  a  function  of  the  Julian  date 
m 

in  the  Smithsonian  Meteorological  Tables  (List,  1958).'  The 
following  table,  after  Manabe  and  Moller  (1961)  was  employed 
to  derive  the  daily  mean  cos  0  and  duration  of  insolation 
as  a  function  of  latitude  and  time  of  the  year.   Computa- 
tions of  F  were  based  upon  mean  mid-seasonal  values  of 

Spring  -  April  15th 

Summer  -  July  15th 
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TABLE  VII.   The  seasonal  and  latitudinal  distributions  of 
duration  of  insolation  and  daily  mean  values  of  cos  6. 


Lat  . 

Fractional  Length 

of  Daytime 
Apr      Jul     Oct 

Jan 

Weight 
cos 

Apr 

ed  Mean 
6 
Jul 

Oct 

Jan 

5 

.508 

.517 

.500 

.496 

.625 

.587 

.614 

.591 

15 

.521 

.537 

.492 

.471 

.618 

.601 

.579 

.549 

25 

.533 

.562 

.483 

.450 

.599 

.593 

.524 

.474 

35 

.546 

.596 

.471 

.421 

.558 

.567 

.458 

.393 

45 

.562 

.637 

.454 

.362 

.501 

.521 

.379 

.317 

55 

.596 

.708 

.437 

.321 

.423 

.453 

.282 

.203 

65 

.629 

.837 

.404 

.208 

.345 

.369 

.176 

.106 

75 

.750 

1.000 

.329 

- 

.241 

.311 

.071 

- 

85 

1.000 

1.000 

_ 

_ 

.168 

.318 

_ 

_ 
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TABLE  VIII.   Incoming  solar  radiation  on  a  horizontal 
surface  at  the  top  of  the  tropopause  (ly  day--*-)  . 


Latitude 

Spring 

S  u  mm  e  r 

Fall 

Winter 

0-10 

884 

867 

844 

784 

10-20 

896 

922 

783 

692 

20-30 

889 

952 

696 

571 

30-40 

848 

965 

593 

443 

40-50 

784 

948 

473 

306 

50-60 

702 

916 

339 

174 

60-70 

604 

882 

195 

59 

70-80 

503 

888 

64 

0 

80-90 

468 

908 

0 

0 
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Fall      -  October  15th 
Winter    -  January  15th 
Joseph's  method  (1966)  was  followed  in  partitioning  the 
effective  solar  insolation  into  a  part  subject  to  water- 
vapor  absorption  but  not  to  Rayleigh  scattering,  and  the 
remainder  subject  to  only  Rayleigh  scattering.   The  presence 
of  clouds  introduced  cloud  reflectivities  in  both  portions. 
Thus,  the  solar  insolation  subject  to  tropospheric  absorp- 
tion comprises  wavelengths  :>  0.9  ym  and  was  obtained  from 


FA  =  0.349F 
A 


(36) 


In  the  wavelengths  <  0.9  ym,  water  vapor  absorption  is 
negligible  and  the  solar  insolation  subject  to  Rayleigh 
scattering  was  obtained  from 


F   =  0.651F 


(37) 


B.   DISPOSITION  OF  FLUX  SUBJECT  TO  SCATTERING 

In  treating  the  fraction  of  the  solar  insolation  subject 
to  scattering,  Joseph  (1966),  found  that  the  Rayleigh 
scattering  values  at  sea  level  in  clear  skies  (after  Coulson, 
1959)  could  be  fitted  by  least  squares  to  the  form 


a   =  0.085  +  0.25074  log  (—  sec  6) 


(38) 


where  it  is  the  observed  surface  pressure.   For  a  full-day 


the  value  of  sec  8  in  (38)  was  treated  as  sec  6,  the  time- 
mean  sec  6 ,  which  was  taken  as  the  inverse  of  the  mean 


cos  6  values  from  Table  VII. 
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Another  reflective  parameter  is  the  surface  albedo,  a  . 

g 

After  numerical  experimentation  designed  to  test  for  a 

global  radiative  heat  balance,  a  constant  average  surface 

albedo  of  a   =  0.14  was  deduced.   It  should  be  noted  that 
g 

this  value  of  a   turned  out  to  be  equal  to  that  deduced  by 

g 

Budyko  (1955).   A  discussion  of  the  validity  of  a  constant 
surface  albedo  will  be  reserved  for  the  surface  balance 
section . 

Considering  the  possibility  of  multiple  reflections  be- 
tween earth  and  atmosphere,  each  of  which  diverts  downward 

the  fraction  a   of  the  earth  surface  reflectance  a  ,  the 
s  g 

clear  sky  insolation  at  earth  after  scattering  is  given  by 
Ig0  =  Fs(l-as)[l+asag+— (asag)n+  ---]  (l-ag)     (39) 

or 

For  the  existence  of  a  cloud  cover,  the  Rayleigh  re- 
flection coefficient  of  (38)  was  applied  at  the  cloud  top 
a  =  0.4,  rather  than  (7  =  1.   At  the  cloud  top,  a  unique 

cloud  reflectivity  fraction  R   of  0.5  was  chosen  everywhere 

J  c 

This  reflection  is  basically  a  Mie-type,  diffuse  back- 
scattering  caused  by  water  droplets  in  the  clouds. 

Thus,  at  the  cloud  top,  the  entering  solar-beam  insola- 
tion is  given  by  n  , 

F„[l-a.(~  sec  8)](1-R  ) 
F     =   s     S   Po c  ,41* 

FSC4  l-a_R  C   ; 

S  c 

This  downward  scattered  beam  is  forced  into  a  diffuse 

angular  array  by  the  Mie  scattering  process.   In  order  to 
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account  for  further  attenuation  by  Rayleigh  scattering  when 
the  ray  paths  are  diffuse,  a  mean  slant  path  angle  6   such 
that  sec  8   =  5/3  was  judged  more  appropriate  than  the  actual 
sec  G  (Katayama,  1966). 

The  Rayleigh  scattering  coefficient  beneath  the  cloud 
top  (i.e.  in  layer  (4,10))  was  found  to  be 


aSD  "  aSDir,  "  aSD/,  "  °-25074  loS  2-5 

10  4 


(42) 


Thus,  a    below  the  cloud  top  has  a  constant  value  of  0.0998. 
In  this  connection,  it  may  be  noted  that  Manabe  and  Strick- 


ler  (1964)  employed  a  constant  value  of  a   =  0.07  S  cos  9 
loss  to  space  by  Rayleigh  scattering  regardless  of  the  pres- 
ence of  clouds. 

Multiple  reflections  take  place  below  the  cloud  top. 
Martin  (1972)  has  shown  that  upon  simplification,  the  inso- 
lation at  earth  after  cloudy-sky  scattering  is 


Icr=  Fcr. (1-a   ) (l-a  )  + 
SC    SC4     SD      g 


FSC4[aSD+ag(1-aSD)K1-C'sD)<1-0g)Rc 


[l-Rc(l-c.SD)   a  ] 


(43) 


In  the  above  accounting,  certain  small  atmospheric  re- 
flection terms  occurring  in  the  direction  cloud  to  earth 
(or  reverse)  have  been  ignored  as  small  compared  to  the  air 
transmission  factor  (1  -  a   ). 

Finally,  the  weighted  mean  earth-surface  absorption  of 
the  scattered  solar  beam  becomes 


Is(10)  =  (1-CL)IS0  +  ISC(CL) 
where  I    is  given  in  (40)  and  I    in  (43). 


(44) 
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C.   TROPOSPHERIC  ABSORPTION  OF  SOLAR  RADIATION 

In  considering  the  absorption  of  the  F.  component  of 
solar  insolation,  the  Manabe-Moller  solar-absorptivity 
function 

a  =  0.0946  (u2  2R  sec  6)0,303  (45) 

was  used.   The  form  of  Equation  45  indicates  that  the  ab- 
sorption function  a  is  applied  over  the  pressure-scaled 
water  vapor  mass  in  the  layer  (2,2K)  along  the  zenith 
path  of  slant  angle  0.   Restricting  absorption  to  the  ratio 
0.349F  requires  that  the  multiplicative  constant  0.0946  of 
(45)  be  adjusted  to  (0.0946/0.349)  =  0.271. 

The  absorbed  insolation  in  the  layer  (2,6)  is  then 
A(2,6)  =  0.271[u(2,6)sec  6]0,303  F  (46) 

Similarly,  the  absorbed  insolation  in  the  layer  (2,10)  is 

A(2,10)=  0.271[u(2,10)sec  9]°,303FA  (47) 

Therefore,  the  absorbed  insolation  in  the  layer  (6,10)  is 
given  by 

A(6,10)=  A(2,10)  -  A(2,6)  (48) 

Since  the  amount  of  solar  beam  absorbed  in  the  layer  (2,10) 
is  available  by  (47),  the  solar  insolation  in  the  F  -parti- 
tion received  at  earth  is  given  by  F  -  A(2,10).   This  por- 
tion is  reduced  by  ground  albedo  a   so  that  the 

g 

earth-absorbed  insolation  resulting  from  F   is 

IAn(10)  =  FA{l-0.271[u(2,10)sec  9 ] ° * 303 } (1-a  )   (49) 
AO  A  g 
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A  modification  is  made  for  cloudy  skies.   The  absorption 
in  the  layer  (2,4)  was  computed  according  to  the  Manabe- 
Moller  absorptivity  function  [Equation  (46)].   The  solar 
beam  at  level  0.4  is  given  by 

FA(4)  =  FA  -  A(2,4)  (50) 

Of  the  F  -energy  available  at  level  O    =  0.4,  the  fraction 
(1-R  )  is  transmitted  through  the  cloud  top.   Thus 

FA(4)  -  [FA-A(2,4)](l-Rc)  (51) 

Within  the  cloud  layer,  the  water  vapor  path  below  a  =  0.4 
is  subject  to  the  mean  diffuse-path  augmentation  factor  5/3 
(after  Katayama,  1966). 

Inside  the  cloud  layer,  another  variable  has  been  intro- 
duced.  The  value  is  u   or  equivalent  water  vapor  mass  for 

c  r 

droplet  absorption  in  clouds.   The  values  of  the  absorption 
coefficient  a     due  solely  to  cloud  particles  have  been 

Lj  Lt  U 

taken  after  theoretical  computations  of  Korb,  Michalowski 
and  Moller  (1956).   These  values  are  listed  after  Manabe 
and  Strickler  (1964,  p.  367)  in  Table  IX,  below. 


TABLE  IX.   Absorptivities  (a    )  of  cloud  droplets. 

L»  J-j  JJ 

Latitude  Band  a 

0-39N,  low-to-middle  clouds  0.035 

40-64N,  low-to-middle  clouds  0.020 

65-90N,  ice-crystal  clouds  0.005 


Since  only  a  single  cloud  layer  was  used  in  the  radia- 
tive model,  the  larger  absorptivities  associated  with  the 
larger  water  content  was  identified  with  the  0-39N  latitude 
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band  of  Table  IX.   The  values  of  a     are  due  solely  to  ab- 

I*  J_t  D 

sorption  by  either  liquid  water  or  the  solid  ice  content  of 
clouds.   The  values  of  a.     are  considered  mean  values  for 

\j  Li  U 


the  latitude  bands  shown  in  Table  IX.   Moreover,  the  a 


CLD 


values  of  Table  IX  are  averaged  for  the  entire  solar  beam, 
and  not  merely  for  the  F   portion  of  the  insolation.   Hence, 
the  value  of  a     has  been  multiplied  by  (1.0/0.349)  for 

KjXj  is 

application  to  the  F  -partition.   In  this  spectral  region, 
the  effective  cloud-water  absorp t ivit ie s  become 


a   =  2.8653  a 


CLD 


(52) 


or 

a*  =  {|(u(4,8)+uc)]°*303-[|u(4,8)]0-303}0.271     (53) 

* 
Equation  (53)  may  be  solved  for  u   as  a  function  of  a   and 


u(4,8) 


u 


u(4,8) 


=  {[1  + 


a   .       0.6   0.303,  3.003_1   , 


0.271vu(4,8) 


(54) 


Values  of  u   were  computed  as  a  function  of  u(4,8)  for  the 
latitudes  of  Table  IX  and  are  shown  in  Appendix  Table  B. 
With  u   known,  the  following  calculations  are  made 


A(4,6)=  FA(4){0.271[|  u<4,6)+  y^]0'303} 
A(6,10)=FA(4){0.271[|(u(4,10)+uc)]°,303}-A(4,6) 


(55) 
(56) 


where  F  (4)  is  evaluated  upon  entering  the  cloud  top. 
Finally,  the  insolation  reaching  earth  after  absorption  in 
cloudy  skies  is 


IAC(10)  =  FA(4)  -  A(4,6)  -  A(6,10) 


(57) 
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With  the  surface  albedo  effect  (1-a  )  included,  the  earth- 

O 

absorbed  part  of  I   (10)  becomes 

IAC(10)  =  i;c(10)  [(l-og)/(l-agRc)]  (58) 


D.  NET  INSOLATION  ABSORBED  AT  EARTH 

The  weighted  sum  for  the  insolation  at  earth  after  ab- 
sorption is  given  by 

IAC(10)  =  CL[IAC  of  (58)]  +  (1-CL)[IA0  of  (49)]  (59) 

Thus  the  total  insolation  absorbed  at  earth  is  given  by 

Q1Q  =  IA(10)  +  V10)  (60) 

The  resulting  values  of  Q  „  are  listed  in  Table  X. 

E.  COMPARISON  OF  MERIDIONALLY  AVERAGED  RADIATION 
QUANTITIES 

In  order  to  compare  the  various  annual  outputs  of  this 
radiation  model  with  the  results  of  other  researchers,  the 
radiational  quantities  derived  by  the  NPS  model  were  time- 
averaged  and  space-averaged  over  the  Northern  Hemisphere. 
The  resulting  global  mean  annual  quantities  are  compared 
with  the  results  by  Katayama  (1967),  Budyko  (1956)  and 
London  (1957).   The  mean  values  are  presented  as  a  percent- 
age of  the  solar  radiation  at  the  top  of  the  troposphere. 

Mer idionally-averaged  values  were  generated  using  the 
formulation 

*2       2 

2tt    /         a      R    cos    cj>    d<j> 

R   =   ± (61) 

*2       2 
2tt    /         a       cos    4>    d(J> 
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TABLE  X.   Seasonal  distribution  of  insolation  absorbed  at 
earth  (Q10)  in  ly  day~l. 


Latitude 

Spring 

Summer 

Fall 

Winter 

0-10 

451 

430 

424 

409 

10-20 

484 

471 

404 

385 

20-30 

4S7 

510 

374 

316 

30-40 

450 

525 

315 

234 

40-50 

405 

484 

241 

157 

50-60 

361 

451 

166 

88 

60-70 

317 

430 

95 

30 

70-80 

264 

430 

31 

0 

80-90 

248 

454 

0 

0 
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where 

R  =  annual  mean  radiative  quantity  at  latitude  (J) 

a  =  radius  of  earth 

$  =  center  of  the  latitude  band  under  consideration 

d<j)  =  latitude  increment  in  degrees  between  successive 
bands.   This  increment  is  a  constant  of  10°  for 
this  study. 

The  formulation  in  Equation  (61)  may  be  reduced  to  a 

series  summation  using  cosine  weighting  factors  as  indicated 

in 


R  cos  $ 


R  = 


n  =  l 


9 
n  =  l 


(62) 


cos  <p 


where  $      =  5°,  4>?  =  15°,  etc. 

When  the  results  of  this  study  are  compared  with  those 

of  Katayama,  Budyko,  and  London  in  Table  XI,  it  is  found 

that  they  are  in  rather  good  agreement  as  a  whole.   However, 

there  are  some  minor  differences.   First,  the  older  model 

estimates  of  solar  radiation  absorbed  in  the  atmosphere  are 

generally  3  to  5  units  higher  than  determined  by  this  study. 

Also,  while  the  average  global  albedo,  A  (<£),  determined  in 

this  study  is  typical  of  usual  values  quoted,  namely  34%, 

the  meridional  distribution  of  A  (d))  was  more  nearly  uniform 

P 

with  respect  to  cf>  than  would  have  resulted  if  the  means  had 
been  based  upon  hourly  computations  for  the  months  involved. 
This  result  leads  to  a  lower  than  normal  value  of  solar 
radiation  absorbed  at  the  surface  in  low  latitudes  and  a 
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TABLE  XI.   Comparison  of  me r idionally  averaged  values  of 
radiation  quantities.   The  values  are  presented  as  per- 
centages of  the  global  mean  insolation  at  the  top  of  the 
troposphere . 

Katayama   Budyko   London   Plante 


Solar  radiation  at 
top  of  troposphere 


100 


100 


100 


100 


Solar  radiation  absorbed 

in  atmosphere  18 

Global  albedo  37 


16 
40 


14 
35 


13 

34 


Solar  radiation  absorbed 
at  surface 


45 


44 


47 


52 


Outgoing  long-wave 
radiation  at  200  mb 


63 


60 


62 


64 


Outgoing  long-wave 

radiation  at  surface  18 

Sensible  and  latent 

heat  flux  (see  Section  IX)    29 


17 


27 


18 
29 


18 


31 
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larger  than  typical  value  at  the   surface  in  high  latitudes. 
The  reason  for  this  is  basically  threefold:   First,  atmos- 
pheric absorption  by  dust  was  not  considered  in  the  NPS 
model,  since  it  is  essentially  a  predictive  model.   However, 
London  (1957)  applied  the  results  of  Houghton  (1954)  who 
used  atmospheric  transmissivity  due  to  dust  as  given  by 

td  -  (0.95)m  (63) 


where  m  =  sec  0,  represents  the  air  mass  traversed  by  the 
solar  beam.  London's  depletion  by  atmospheric  dust  would 
then  be  given  by 

D  =  (1  -  TD)  (64) 

of  which  25%  is  considered  subject  to  absorption.   Even 
so,  the  contribution  of  dust  absorption  is  small  and  can 
generally  be  considered  of  the  order  of  1%,  which  is 
approximately  the  difference  between  the  NPS  absorption 
figure  and  that  of  London's. 

Another  contribution  to  absorption  not  previously  con- 
sidered here  is  secondary  absorption  of  the  reflected  por- 
tion of  the  solar  beam.   Since  the  reflected  portion  of  the 
solar  beam  above  clouds  is  subject  to  the  same  percentage 
absorption  as  the  incoming  beam,  the  error  in  absorption 
may  be  as  high  as  1%  of  the  incoming  solar  radiation.  Neither 
London's  study  nor  this  study  considered  secondary  absorp- 
tion of  the  reflected  beam  and  consequently  the  resulting 
atmospheric  absorptions  are  closely  comparable.   Katayama 
considered  secondary  absorption  following  a  first  reflection 
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of  the  solar  beam  by  a  cloud  or  by  ground  or  both.   As  a 
result  Katayama's  solar  absorption  in  air  was  the  largest 
of  the  four  cases  considered. 

The  third  consideration  is  Katayama's  cloud  cover.  In 
the  Katayama  radiative  study,  a  larger  climat ological  cloud 
cover  was  considered  than  in  the  present  study,  or  that  of 
London  or  Budyko .   The  cloud  amounts  were  generally  charac- 
teristic of  coverages  over  the  oceans  (after  McDonald,  1938) 
rather  than  a  combined  ocean-continent  regime. 

It  is  not  exactly  clear  how  Katayama  allowed  McDonald's 
overlapping  cloud-layers  to  be  combined,  although  he  con- 
strained the  total  cloud-coverage  to  agree  with  that  of  Sir 
Napier  Shaw's  (1936).   The  significant  feature  of  the 
greater  cloud  coverages  is  that  they  afford  not  only  greater 
solar  absorptivity  but  also  a  greater  global  albedo  as 
postulated  by  Katayama's  empirical  formulas 


Ac  =  A^  [l-exp(-b1  h)] 


Rc  =  R^  [l-exp(-b2  h)] 


(65) 
(66) 


where  h  is  the  cloud  thickness  and  b,,  b»,  A  ,  R   are 
characteristic  values  for  a  given  cloud  type.   Thus  Kata- 
yama finds  a  greater  cloud  absorptivity  than  was  found  by 
London  or  Plante  (Table  XI)  as  well  as  a  larger  global 
albedo  because  of  the  larger  cloud  coverages. 

F.   GLOBAL  ALBEDO  CONSIDERATIONS 

In  the  process  of  deriving  the  transmission  of  the  two 
solar  components  F   and  F  ,  the  following  dispositional 
terms  were  discussed: 
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(i)     the  combined  insolation  rate  at  the  tropopause 
denoted  by  F  of  Equation  (35). 

(ii)    the  absorption  rate  of  insolation  in  the  atmos- 
phere 

(iii)   the  absorption  rate  of  insolation,  Qin)  at  earth 

With  these  terms,  the  reflected  insolational  rate  from 

planet  earth,  as  a  function  of  latitude  and  season  may  be 

computed.   This  rate  may  be  denoted  by  A  F,  where  A   is  the 

P  P 

global  albedo.   The  term  A  F  must  be  identified  as  the 


difference  of 


F  -  A(2,10)  -  Q1Q  =  ApF 


(67) 


It  should  be  noted  that  three  significant  reflective 
parameters  have  affected  the  value  of  the  left  side  of  (67) 
namely  the  sky  albedo,  the  cloud  albedo  (R   =  0.5),  and  the 
earth-surface  albedo,  a  .   For  simplicity,  the  latter  two 
were  held  constant  at  reasonably  well  accepted  figures  in- 
dependent of  latitude.   While  direct  computations  of  the 
right  side  of  (67)  were  made  at  each  step  of  computing  (ii) 
and  (iii),  the  resulting  global  albedo  values  must  also  be 
consistent  with  the  values  inferred  from  the  left  side  of 
(67). 

A  note  of  caution  in  regard  to  comparing  global  albedo 
as  shown  in  Table  XI.   A  correct  mean  global  albedo  may 
distort  the  atmospheric  heat  balance  if  not  distributed 
correctly  with  respect  to  latitude.   For  example,  if  a 
meridional  bias  exists  in  the  computed  global  albedo,  A  , 
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with  too  large  a  value  at  low  latitudes  and  too  small  a 
value  at  high  latitudes,  the  result  will  influence  the  pole- 
ward heat  transport  required  for  balance,  after  independent 
computation  of  latent  and  sensible  heat  transfers  across 
the  earth-air  interface  (see  Section  IX) .   This  type  of 
bias  in  global  albedo  appears  to  have  developed  in  this 
study  presumably  by  oversimplifying  the  ray-paths  in  clouds 
Choosing  a  single  mean  path  through  clouds  for  a  given  re- 
presentative day,  and  time  of  year  seems  to  have  overesti- 
mated the  global  albedo  in  this  study  in  lower  latitudes 
while  underestimating  it  at  high  latitudes. 
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V.   RADIATIVE  TEMPERATURE  CHANGES 

Long-wave  radiation,  with  a  normal  temperature  lapse, 
will  lead  to  a  cooling  of  the  atmosphere  as  seen  as  seasonal 
distributions  in  Tables  v(a-d).   In  certain  locations  and 
seasons,  some  of  the  atmospheric  layers  considered  were 
subject  to  temperature  inversions  which  led  to  local  heating 
rather  than  the  usual  cooling.   Another  heating  situation 
examined  was  at  an  overcast  cloud  base.   The  terminology, 
'cloud  base  warming'  or  'green  house  effect'  may  be  used 
to  describe  this  phenomenon.   Since  the  net  flux  within  a 
cloud  layer  is  zero,  a  negative  net  flux  increment  existed 
in  the  cloud  to  sub-cloud  layer  (8,10).   With  the  sign 
convention  adopted  here,  this  represents  a  base-warming 
effect.   In  the  cloud  top   region,  there  is  a  substantial 
net  flux  increment  F2   -  F,   which  leads  to  considerable 
cloud-top  cooling.   Thus,  in  infrared  radiative  temperature 
changes  the  role  of  clouds  is  clearly  to  increase  the  cooling 
rate  at  the  cloud  tops  and  to  decrease  the  cooling  rate 
at  the  cloud  bases.   This  process  can  be  considered  as  a 
decrease  of  the  atmospheric  stability  within  the  cloud  layer. 

The  mean  annual  infrared  temperature  changes  in  the 
troposphere  are  plotted  in  Figure  7  as  a  function  of  latitude 
The  dependency  of  these  cooling  rates  upon  the  climato- 
logical  cloud  cover  at  these  same  latitudes  is  evident  by 
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inspecting  Figure  8.   A  linear  regression  was  performed  re- 
lating the  mean  long-wave  temperature  changes  of  each  tropo- 
sphere layer  against  the  total  cloud  cover  used  in  the 
computation  at  latitude  <J> .   Each  latitude  band  for  each 
season  was  considered  separately,  giving  a  total  of  36  cases 
For  the  layer  (6,10)  the  correlation  coefficient  was  -0.95 
while  for  the  layer  (2,6)  the  correlation  was  +0.99.   Graphs 
showing  the  comparison  of  temperature  changes  found  with 
clear  and  with  average  cloudy  skies  for  the  layer  (2,6)  and 
layer   (6,10)  are  shown  in  Figures  9  and  10,  and  clearly 
demonstrate  the  phenomena  of  cloud-top  cooling  and  cloud- 
base  warming. 

In  conditions  of  normal  cloudiness,  the  infrared  heat 
losses  are  partially  offset  by  atmospheric  absorption  of 
solar  radiation.   The  resultant  temperature  change  in  the 
layer  (6,10)  is 

it  ~(F**-  Fin*)  +  A(6>in> 

If  <6'10>  "     0.4,(1000/g)  c (68> 


The  temperature  change  in  the  layer  (2,6)  is  given  by 

"(F2*~  F6*}  +  A(2'6) 


fl  <6>10>  " 


0.4TT(1000/g)  c. 


(69) 


In  most  cases,  this  absorption  of  solar  insolation  is  smal- 
ler than  the  cooling  that  results  from  long  wave  radiation 
and  acts  only  to  modify  the  pattern  of  radiative  heat  loss, 
as  shown  in  Figure  11.   The  dashed  curves  represent  the 
annual  mean  rate  of  temperature  change  due  to  long  wave 
radiative  cooling  for  the  layers  (6,10)  and  (2,6).   The 
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solid  curves  represent  the  result  of  net  radiative  cooling 
(combined  effects  of  long  and  short  wave  radiation)  in  the 
same  layers.   Atmospheric  absorption  clearly  only  modified 
the  magnitude  of  radiative  cooling  while  the  shape  of  the 
curves  was  unchanged.   By  comparison  of  Figures  8  and  11, 
high  correlations  between  net  radiative  temperature  changes 
and  total  cloud  cover  may  once  again  be  noted. 
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VI .   RADIATIVE  BALANCE  AT  THE  TROFOPAUSE 

A  radiative  budget  was  calculated  at  the  tropopause  and 
compared  to  the  results  of  Houghton  (1954).   Houghton  made 
a  comparison  of  the  net  incoming  solar  radiation  at  the 
tropopause  with  the  outgoing  long-wave  flux.   He  determined 
that  insofar  as  the  earth-atmosphere  combined  system  is  con- 
cerned, a  column  at  about  35  N  latitude  experiences  a  mean 
radiative  energy  balance.   Corresponding  latitude  belts  of 
the  system  located  equatorward  of  latitude  35N  have  a  mean 
annual  radiative  energy  surplus.   Those  latitude  bands  pole- 
ward  of  35  N  latitude  experience  a  mean  annual  radiative 
energy  deficit.   These  mean  annual  radiative  surpluses  and 
deficits  impose  requirements  for  horizontal  heat  transfer 
rates  across  latitude  walls  in  both  air  and  ocean. 

The  net  insolation  at  the  tropopause  (Q  )  is  given  by 


%  "   «I  "  QR 


(70) 


where  Q   is  the  incident  insolation  at  the  tropopause  after 


a  4%  reduction  due  to  attenuation  in  the  stratosphere.   Q 

is  the  reflected  insolation  returned  to  space.   The  ratio 

of  Q  /[1.04  Q  ]  is  the  measure  of  the  earth  atmosphere 
R  I 


R 


planetary  albedo.   Q   is  given  by 


«R  =  1l  -  QA  "  Q10 


(71) 


where  Q   is  the  short  wave  radiation  absorbed  in  the  atmos- 

A. 

phere  and  Q  _  is  the  total  short  wave  radiation  absorbed 
at  earth. 
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The  resultant  Q   is  a  function  of  the  Rayleigh  back 
scattering,  cloud  reflectivity  and  surface  albedo.   The 
factors  determining  Rayleigh  scattering  were  discussed 
in  Section  IV. 

Treatment  of  cloud  reflectivity  was  limited  here  by 
the  use  of  a  single  cloud  layer.   A  global  average  of 
about  one-fourth  of  the  incident  solar  radiation  is  re- 
flected to  space  by  clouds.   In  this  model  a  single  cloud 
albedo  R   =  0.5  was  used.   Sellers  (1965)  employed  a  range 
of  albedos  for  differing  cloud  types.   A  similar  treatment 
using  various  cloud  types  with  various  albedos  was  adopted 
by  London.   On  a  mean  annual  basis,  0.5  appears  to  be  a 
good  approximation  for  the  hemispheric  cloud  reflectivity, 
but  on  an  individual  gridpoint  basis  the  cloud  albedos 
adopted  in  this  study  may  be  an  oversimplification.   Future 
testing  in  this  area  is  envisioned  for  the  NPS  predictive 
model . 

A  single  constant  surface  albedo  of  0.14  was  chosen 
after  experimentation  with  the  tropospheric  energy  balance. 
Budyko  (1955)  estimated  the  mean  earth  surface  albedo  as 
0.14  while  Katayama  (1967)  found  a  value  of  0.13.   But  the 
albedo  varies  with  the  type  of  surface.   The  short-wave 
albedo  is  highest  over  fresh  snow  (  90%)  and  is  lowest  over 
water  and  forests  (approximately  5%) .   Since  there  are 
seasonal  as  well  as  annual  variations  in  the  extent  of 
vegetation  and  snow  cover,  an  attempt  at  strict  evaluation 
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of  ground  albedo  for  these  variations  is  at  best  difficult. 
Thus,  for  simplicity,  a  single  mean  ground  albedo  was  em- 
ployed . 

The  net  outgoing  long-wave  radiation  at  the  tropopause 
was  given  by  F„   since  the  0.2  sigma  level  was  adopted  as 
the  mean  height  of  the  tropopause. 

Therefore,  the  net  radiative  flux  across  the  tropopause 
is  given  by 


%    '    "o  -  F2 


(72) 


Results  for  seasonal  variations  in  the  net  flux  across  the 
tropopause  are  listed  in  Tables  Xll(a-d).  The  mean  annual 
radiative  input  by  latitude  bands  is  listed  in  Table  XII(e) 
along  with  the  cosine-latitude  weighting  factors.  The  re- 
sulting mean  Northern  Hemisphere  radiative  budget  is  zero 
as  anticipated,  since  it  is  assumed  that  the  hemisphere  is 
in  annual  radiative  balance. 

During  the  summer,  the  tropopause  has  a  net  radiative 
surplus  throughout  all  latitude  bands.   As  seen  in  Table 
XII(b),  the  incident  insolation  is  largest  during  this  sea- 
son and  far  exceeds  the  small  increase  in  net  long  wave 
radiation  due  to  warmer  summer  temperatures.   The  opposite 
case  is  true  for  winter.   Low  insolation  rates  produced  a 
negative  radiative  budget  at  all  latitude  bands  except  the 
0-10°  band.   This  particular  latitude  band  has  a  positive 
radiative  net  surplus  during  the  entire  year.   Spring  and 
fall  show  transitions  between  the  extremes  of  summer  and 
winter . 
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Effects  of  variations  in  cloud  cover  on  the  computed 
net  radiative  input  at  the  tropopause  can  be  seen  in  Table 
XIII  which  lists  spatial  mean  averages  of  the  radiative 
quantities  across  the  tropopause.   It  is  seen  that  clouds 
exert  an  overall  cooling  effect  in  this  model. 


TABLE  XIII.   Spatial  mean  annual  net  radiative  input  at 
the  tropopause  (Northern  Hemisphere). 


CL  =  0    CL  =  1.0    Average  CL 


694 
130 
564 
528 


694 
360 
334 
373 


694 
244 
450 
450 


Net  Tropopause 
Radiative  Input 


36 


-39 
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VII  .   GLOBAL  ALBEDO 

The  earth's  global  albedo  results  from  the  upward  re- 
flection by  the  ground,  by  clouds  and  by  scattering  by 
atmospheric  molecules  of  solar  radiation.   Variations  in 
computed  planetary  albedo  are  strongly  dependent  on  total 
cloud  cover.   The  distribution  with  latitude  and  season 
of  global  albedo  is  listed  in  Table  XIV.   There  is  a  pro- 
nounced minimum  in  the  subtropics,  which  is  a  region  of 
minimum  cloudiness.   The  minimum  shifts  equatorward  during 
the  winter  and  poleward  during  the  summer,  following  move- 
ment of  the  ITCZ.   The  increased  albedo  in  high  latitudes 
during  all  seasons  results  from  an  increase  in  cloudiness 
with  latitude  up  to  about  60  to  70N. 

Total  global  albedo  when  averaged  over  all  latitudes 
and  seasons  is  about  34  percent.   The  total  global  albedo 
shows  little  variation  with  season  and  agrees  closely  with 
London's  results.   Figure  12  shows  a  comparison  between 
mean  annual  planetary  albedo  with  latitude  computed  by 
London  and  the  results  of  this  study.   Even  though  the 
spatially  averaged  annual  global  albedo  of  34%  agrees  close- 
ly with  London,  there  is  disagreement  in  its  distribution 
with  latitude,  a  fact  referred  to  in  Section  IV. 

In  tropical  latitudes,  London's  global  albedo  was  less 
than  that  computed  in  this  study,  while  in  the  northern 
latitudes,  his  values  exceeded  those  computed  here.   Since 


76 


\ 


\ 


\ 


\ 


o 


< 


4-» 


CO 

O 
•H 
4-1 

-Q 
•H 
M 
4-» 
W 
•H 

O 

CD 
^3 
rH 

CO 

H 
ctf 

,Q 
O 

rH 

60 


ctj 

c 

CO 

ctj 

0) 

0 

•4-1 

o 

o  • 
w   <y 

•H  'V 
U    3 

Ctj     4J 

G.  -H 
0  *J 
O  Ctj 
U  H 


CN 


U 

bO 
•H 
p4 


O 

m 


o 


o 

CO 


(  V)  opaqiV  I^qoTO  luaojsj 


77 


TABLE  XIV.   Seasonal  distribution  of  planetary  albedo  with 
latitude  (in  percent  of  QT). 


Latitude 

Spring 

Summer 

Fall 

Winter 

0-10 

35 

37 

36 

34 

10-20 

32 

35 

34 

30 

20-30 

32 

32 

32 

31 

30-40 

36 

32 

34 

36 

40-50 

39 

37 

38 

39 

50-60 

40 

40 

41 

41 

60-70 

39 

40 

42 

40 

70-80 

40 

41 

41 

- 

80-90 

40 

40 

- 

— 
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identical  total  cloud  cover  amounts  were  used  in  both 

studies,  the  treatment  of  atmospheric  and  cloud  albedos 

seem  to  be  the  major  cause  for  this  difference.   The 

present  study  used  a  constant  cloud  albedo  of  R   =  0.5 

c 

while  London  used  varying  cloud  albedos  associated  with 
constituent  cloud  types  after  Haurwitz  (1948).   London 
thus  used  a  multi-layer  cloud  model  while  this  was  not 
done  here.   Furthermore,  the  device  of  replacing  the 
diffuse  path  length  in  clouds  by  an  equivalent  path  5/3 
longer  seems  to  give  too  short  a  path  subject  to  cloud- 
reflecting  processes  at  high  latitudes  and  too  long  a 
path  at  low  latitudes.   This  is  important  if  cloud-reflec- 
tance is  to  be  treated  by  Mie-scat ter ing  theory.   Varia- 
tions of  surface  reflectance  with  solar  zenith  angle  also 
plays  a  secondary  role  in  explaining  the  variations  of 

A   with  latitude. 
P 
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VIII.   ATMOSPHERIC  BUDGET 

A  simplified  annual  atmospheric  heat  budget  for  the 
Northern  Hemisphere,  along  with  seasonal  distributions  of 
radiative  surpluses  and  deficits  was  constructed.   The 
stratosphere  was  assumed  to  be  in  radiative  equilibrium. 

The  atmosphere  radiation  balance  is  shown  schematically 
below : 

«i         Or  f*2 

_1 1 L_ 

Qio         qe     Qs         »io 


1 _J L 


The  net  cooling  contribution  due  to  long  wave  radiation  is 
given  by 

FN  "  F2  *  F10  <73) 

The  effective  insolation  at  the  top  of  the  troposphere 
is  the  total  incident  solar  radiation  minus  the  reflected 
insolation  to  space.   The  total  atmospheric  absorption  is 
therefore 

QA  =  (Qj.  -  <y  -  Q10  (74) 

where  Q    is  the  short  wave  radiation  absorbed  at  earth 
using  a  mean  surface  albedo,  a   =  .14.   A  final  heating 
term  consists  of  the  latent  heat  of  evaporation  and  sen- 
sible heat  transfer  rates  across  the  earth-air  interface 


and  denoted  by  Q  : 

H 


QH  =  QE  +  Qs  (75) 
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QVA  "  ^H  +  «A  -  FN 


The  parameterization  of  Q£  and  Q   is  detailed  in  Section  IX. 
For  each  latitude  band,  the  net  atmospheric  heating  rate  by 
the  three  mechanisms  (insolation  absorption,  long-wave 
absorption,  interface  convection)  is  given  by: 

(76) 

The  resulting  atmospheric  heating  (cooling)  rates  in  ly  day 
are  summarized  in  Tables  XV(a-d)  for  each  latitude  band  and 
each  of  the  four  seasons,  while  Table  XV(e)  gives  the  annual 
mean,  derived  by  averaging  the  four  seasonal  heating  rates, 
as  a  function  of  latitude. 

It  is  to  be  noted  that  QTT  .  has  been  termed  the  "balance 

VA 

requirement"  by  Davis  (1963),  and  represents  the  excess 
amount  to  be  disposed  of  from  a  latitude  belt  by  poleward 
flux-divergence  of  heat.   Table  XV(e)  shows  that  the  global 
annual  average  of  this  parameter  (as  calculated  in  this 
study)  is  zer o . 

Warming  of  the  atmosphere  by  the  absorption  of  solar 
radiation  is  less  than  the  long-wave  cooling  rate.   The  net 
result  is  that  the  atmosphere  is  cooling  radiatively  at  all 
latitudes  at  a  rate  of  about  1.1°C  day   .   The  compensating 
energy  requirement  for  possible  balance  has  been  indicated 
by  Equation  (7 6) as  arising  by  sensible  and  latent  heat  trans- 
fers across  the  earth  surface. 

With  the  addition  of  sensible  and  latent  heat  into  each 
latitude  band,  Table  XI(e)  shows  that  latitude  bands  equa- 
torward  of  40N  have  a  total  energy  surplus  while  latitude 
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bands  poleward  of  AON  exhibit  an  energy  deficit.   While  the 
mechanisms  of  horizontal  transport  are  not  specified  in 
this  studys  the  resultant  globally-weighted  values  of  Q 
of  Table  XV(e)  indicate  a  mean  annual  total  heat  balance 
in  the  atmosphere.   The  poleward  redistribution  required  to 
effect  this  total  heat  balance  is  evidently  brought  about 
by  general  circulation  processes,  which  were  not  considered 
in  this  study.   For  our  purposes  it  suffices  to  show  that 
the  energy-generating  aspects  of  the  model  proposed  here  is 
subject  to  an  atmospheric  balance. 
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IX.   SURFACE  HEAT  TRANSFERS 

The  earth's  surface  heat  balance  was  computed  using 
evaporation  and  sensible  heat  transfers  over  the  oceans 
as  surface  heat  sinks.   Heat  transfer  due  to  evaporation 
and  sensible  heat  are  more  easily  computed  over  ocean 
surfaces  using  climatological  data  available  in  the  Paci- 
fic Ocean  representative  of  the  latitude  belt  as  a  whole. 
In  effect,  London's  atmospheric  soundings  were  used  in 
association  with  data  for  ocean  surface  temperatures  and 
surface  winds  taken  from  U.  S.  Weather  Bureau's  Atlas  of 
Climatic  Charts  of  the  Oceans  (1938).   Continuity  of 
sounding  data  with  the  oceanic  reports  was  maintained  by 
maintaining  the  same  temperature  lapse  rate  in  the  layer 
zero  to  1  km. 

The  surface  radiative  budget  is  shown  schematically 
below : 


Qio    Fio 


%         Qs 


It        ft 

and  may  be  stated  symbolically  by  the  equation: 


Qm  "  <F^   +  Qw  +  Qq>  =  Q 


10 


10 


E    ^S 


VO 


(77) 


Here  Q  is  the  required  oceanic  heat  transport  divergence. 
The  surface  layer  evaporation  heat  transfer  rate  was  calcu- 
lated after  Langlois  and  Kwok(1969)  using: 

%  '    L  p10  CD  VS  ("sg  "  «,>  (78) 
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with  the  constants: 


L  =  590  cal  gram 


-1 


-3         -3 

p_n  =  1.2  x  10    gram  cm 

Also  V   is  the  climatological  ocean  surface  wind  speed  in- 
cremented by  2.2  mps  in  order  to  account  for  the  effect 
of  gustiness  upon  evaporation  [Langlois  and  Kvok  (1969)]. 
A  variable  drag  coefficient  C   proposed  by  Deacon  and  Webb 

(1962)  for  the  case  of  neutral  stability  was  employed: 

-3 


Cn  =  (1.0  +  0.07  V  )  x  10 


(79) 


At  the  ocean  surface,  the  atmosphere  was  assumed  to 

be  saturated  with  respect  to  a  plane  water  surface.   The 

surface  specific  humidity,  q   ,  was  computed  using  the 

sg 

Claus ius-Cl apey ron  equation: 

(AE-BE/TG) 


.622  e 


sg 


TT-e 


(AE-BE/TG) 


(80) 


The  specific  humidity  at  the  top  of  the  constant  flux  sur- 
face layer,  q  ,  is  given  following  the  parameterization  of 
the  planetary  boundary  layer  by  Langlois  and  Kwok,  by  the 
formula : 


kq9  +  (y  Z1Q)   Cp  Vs  qsg 
k  +  <V  Z10>  S  VS 


with: 


k  = 


a  ( 


(81) 


(82) 


1  + 


Z9  "  Z10 


*       * 

In  (82),  k   and  a   are  turbulent  transfer  coefficients 

whose  values  have  been  documented  by  Kaitala  (1972). 
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Ice  free  oceans  exist  equatorward  of  60N  for  all  seasons, 
Therefore,  air  and  sea  temperature  and  wind  data  were  re- 
quired for  only  the  first  six  latitude  bands.   Since  Lon- 
don's atmospheric  data  was  distributed  at  one-kilometer 

intervals,  (Z   -  Z   )  was  taken  as  1  km;  also  T   and  q   are 
y    j.  u  y      y 

temperature  and  specific  humidity  at  the  one-km  height.  The 
first  kilometer  temperature  lapse  rate  was  computed  using 
London's  soundings.   Air  temperature  at  the  ocean  surface 
was  assumed  equal  to  the  ocean  surface  temperature.   The 
lapse  rate  was  then  applied  to  yield  a  new  TQ  for  use  in 

y 

(82).   The  new  specific  humidity  q   at  1  km  was  generated 
using  the  relative  humidity  from  London's  data  and  the 
Clausius  Clapeyron  equation. 

Sensible  heat  transfer  rates  equatorward  of  60N  were  com- 
puted using  Bowen  ratio  (r)  values  from  Sellers  (1965)  and 
the  formulation: 


Qs  =  r   QE 


(83) 


Poleward  of  60N,  combined  evaporation  and  sensible  heating 
rates  were  taken  from  an  assessment  of  many  thousands  of 
ship  (and  also  of  ice  island)  reports  made  by  Vowinckel  and 
Taylor  (1964) . 

Tables  XVI(a-d)  contain  the  seasonal  surface  radiative 
balances,  Q  ,  which  is  the  oceanic  heat  transfer  require- 
ment by  poleward  flux-divergence.  Seasonal  variations  can 
be  compared  against  the  annual  mean  values  listed  in  Table 
XVI(e) . 
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The  physical  meaning  of  a  negative  Q    value  for  a 
latitude  belt  in  the  tropics  (with  positive  values  in 
northerly  belts)  is  that  a  flux-convergence  (divergence) 
of  sensible  heat  due  to  ocean  currents  is  needed  to 
maintain  a  long-term  balance. 
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X.   CONCLUSION 

With  the  use  of  mean  seasonal  atmospheric  data,  a  pro- 
posed radiation-heating  model  for  the  Naval  Postgraduate 
School  primitive-equation  predictive  system  was  evaluated. 
The  model  proved  readily  adaptable  to  computerization. 
The  time  and  space  averaged  radiational  quantities 
showed  rather  good   agreement  as  a  whole  when  compared  with 
the  results  of  other  researchers.   Minor  difference  in  the 
comparisons  can  generally  be  ascribed  to  the  treatment  of 
global  albedo,  and  the  assumption  of  but  a  single  cloud 
layer.  While  the  average  global  albedo  determined  in  this 
study  is  typical  of  usual  values  quoted,  namely  34%,  the 
meridional  distribution  of  global  albedo  was  more  nearly 
uniform  with  respect  to  latitude  than  would  have  resulted 
if  the  mean  value  had  been  based  upon  hourly  computations 
for  the  months  involved.   This  result  led  to  a  lower  than 
normal  value  of  solar  radiation  absorbed  within  the  earth's 
surface  in  low  latitudes  and  a  larger  than  typical  value 
absorbed  in  high  latitudes. 

Of  the  three  heat  budgets  discussed,  results  compatible 
with  Houghton  (1954)  and  Davis  (1963)  respectively,  were 
demonstrated  for  the  earth-atmosphere  system  and  for  the 
magnitude  and  direction  of  the  required  total  heat  flow. 
Similar  tests  were  also  applied  to  the  ocean  surface 
where  the  distribution  of  surface  heat-flux  divergence  was 
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opposite  to  that  expected  (e.g.  Budyko  (1955)),  the  results 
of  the  surface  heat  budget  were  somewhat  biased  by  the 
averaging  involved  in  the  derivation  of  the  global  albedo, 
as  explained  previously. 

Finally,  while  the  radiation  package  has  proven  general- 
ly adaptable  to  a  general  circulation  predictive  model,  fur- 
ther refinements  are  anticipated  in  the  treatment  of  cloud 
cover.   A  two-layer  cloud  model  for  the  treatment  of  middle 
and  low-level  clouds  has  been  suggested.   These  two  cloud 
types  would  have  different  short-wave  flux  reflection  and 
absorption  characteristics  and  would  lead  to  a  more  accurate 
determination  of  global  albedo. 
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APPENDIX    A 
^Evaluation  of   the   final   temperature-dependent   flux   integral 


¥  (v 

Tx)  dB 
and 

for  a  variety  of  final  temperature  values^ (T1 ' 

final  sea 

led  water 

-vapor  values  (ly 

day"1) 

Scaled 

Water 

190° 

195° 

200° 

205° 

210° 

215° 

220° 

Vapor 

0.0001 

26.98 

28.53 

30.12 

31.76 

33.45 

35.19 

36.97 

0.0005 

42.62 

45.36 

48.19 

51.13 

54.16 

57.30 

60.55 

0.001 

50.84 

54.25 

57.80 

61.48 

65.30 

69.26 

73.36 

0.002 

59.89 

64.09 

68.46 

73.01 

77.75 

82.66 

87.78 

0.004 

69.65 

74.74 

80.05 

85.60 

91.38 

97.41 

103.68 

0.006 

75.63 

81.29 

87.20 

93.39 

99.85 

106.59 

113.61 

0.008 

79.97 

86.05 

92.41 

99.08 

106.04 

113.32 

120.91 

0.01 

83.37 

89.79 

96.52 

103.57 

110.94 

118.64 

126.69 

0.015 

89.63 

96.68 

104.09 

111.86 

120.00 

128.53 

137.44 

0.02 

94.09 

101.61 

109.52 

117.83 

126.54 

135.67 

145.22 

0.03 

100.38 

108.57 

117.21 

126.30 

135.83 

145.84 

156.34 

0.04 

104,80 

113.49 

122.65 

132.30 

142.45 

153.11 

164.29 

0.06 

110.94 

120.33 

130.25 

140.71 

151.73 

163.33 

175.51 

0.08 

115.20 

125.09 

135.56 

146.60 

158.25 

170.52 

183.42 

0.1 

118.42 

128.71 

139.60 

151.10 

163.24 

176.04 

189.51 

0.15 

124.08 

135.08 

146.73 

159.07 

172.11 

185.87 

200.37 

0.2 

127.90 

139.40 

151.59 

164.52 

178.19 

192.64 

207.89 

0.3 

132.96 

145.15 

158.10 

171.84 

186.40 

201.81 

218.10 

0.4 

136.30 

148.96 

162.43 

176.74 

191.92 

208.00 

225.01 

0.5 

138.72 

151.74 

165.61 

180.35 

196.00 

212.60 

230.16 

0.6 

140.59 

153.89 

168.07 

183.16 

199.19 

216.20 

234.22 

0.8 

143.31 

157.05 

171.72 

187.34 

203.95 

221.59 

240.30 

1.0 

145.22 

159.29 

174.32 

190.33 

297.38 

225.50 

244.73 

1.5 

148.23 

162.86 

178.50 

195.20 

213.00 

231.96 

252.10 

2.0 

149.99 

164.97 

181.02 

198.17 

216.48 

235.98 

256.73 

2.5 

151.12 

166.36 

182.70 

200.18 

218.85 

238.76 

259.97 

3.0 

151.89 

167.33 

183.90 

201.63 

220.59 

240.81 

262.36 

4.0 

152.82 

168.55 

185.44 

203.55 

222.92 

243.61 

265.68 

5.0 

153.29 

169.23 

186.35 

204.72 

224.38 

245.41 

267.86 

6.0 

153.53 

169.61 

186.90 

205.46 

225.35 

246.63 

269.40 

8.0 

153.59 

169.88 

187.41 

206.24 

226.45 

248.09 

271.23 

10.0 

153.40 

169.81 

187.50 

206.51 

226.94 

248.82 

272.25 

15.0 

152.49 

169.08 

186.98 

206.26 

227.00 

249.26 

273.11 

25.0 

150.35 

167.04 

185.09 

204.58 

225.57 

248.14 

272.38 

30.0 

149.32 

166.02 

184.09 

203.61 

224.66 

247.30 

271.63 
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APPENDIX  B 

Water-vapor  equivalent  u   for  droplet  absorption  in  relation 
to  the  pressure-scaled     water-vapor  absorber  mass  u(4,8) 
(in  the  sounding  layer  (4,8)).   Case  u  (1)  corresponds  to 
low  latitude  water-drop  absorptivity;   u  (2)  to  mid-latitude 
absorptivity;  u  (3)  to  high  latitude  absorptivity. 


u(4,8) 

0.01 
0.02 
0.05 
0.10 
0.15 
0.20 
0.30 
0.50 
0.80 
1.00 
1.20 
1.40 
1.60 
1.80 
2.00 
2.20 
2.40 
2.60 
2.80 
3.00 


u  (1) 

c 

0.14170 
0.18936 
0.28886 
0.40852 
0.50529 
0.56995 
0. 73785 
0.98617 
1.29740 
1.48119 
1.65219 
1.81232 
1.96648 
2 .11306 
2.25403 
2.39016 
2.52206 
2.65022 
2.77499 
2.89675 


u  (2) 
c 

0.05118 

0.97292 

0.12003 

0.17851 

0.22673 

0.26939 

0.34473 

0.47287 

0.63548 

0.73226 

0.82271 

0.90822 

0.98974 

1.06793 

1.14331 

1.21623 

1.28699 

1.35584 

1.42297 

1.48854 


u  (3) 

c 

0.00740 
0.01155 
0.02103 
0.03329 
0.04365 
0.05294 
0.96955 
0.09826 
0.13520 
0.15739 
0.17822 
0.19800 
0.21691 
0.23510 
0.25267 
0.26970 
0.28626 
0.30240 
0.31816 
0.33357 
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APPENDIX    C 


The    distribution    of    relative    humidity    and    total    optical 
mass     (to    the    tropopause)     u(gm    cm-2) (from    London,     1957). 
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APPENDIX    D 

The    distribution    of    average    temperature    and    pressure     (from 
London,     1957) . 
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